Dielectric Susceptibility and Heat Capacity of Ultra-Cold Glasses in Magnetic Field 
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Recent experiments demonstrated unexpected, even intriguing properties of certain glassy mate- 
rials in magnetic field at low temperatures. We have studied the magnetic field dependence of the 
static dielectric susceptibility and the heat capacity of glasses at low temperatures. We present a 
theory in which we consider the coupling of the tunnelling motion to nuclear quadrupoles in order 
to evaluate the static dielectric susceptibility. In the limit of weak magnetic field we find the res- 
onant part of the susceptibility increasing like B^ while for the large magnetic field it behaves as 
1/ B. In the same manner we consider the coupling of the tunnelling motion to nuclear quadrupoles 
and angular momentum of tunnelling particles in order to find the heat capacity. Our results show 
the Schotky peak for the angular momentum part, and B^ dependence for nuclear quadrupoles 
part of heat capacity, respectively. We discuss whether or not this approach can provide a suitable 
explanation for such magnetic properties. 

PACS numbers: 61.43.Fs,77.22.Ch,65.60.-|-a 



I. INTRODUCTION 

At very low temperatures, glasses and other amor- 
phous systems show similar thermal, acoustic and dielec- 
tric properties^, which are in turn very different from 
those of crystaUine solids. Below IK, the specific heat 
of dielectric glasses is much larger and the ther- 
mal conductivity k orders of magnitude lower than the 
corresponding values found in their crystalline counter- 
parts, depends approximately linearly and k almost 
quadratically on temperature?, respectively. This is in 
clear contrast to the cubic dependence observed in crys- 
tals for both properties, well understood in terms of De- 
bye's theory of lattice vibrations. Above IK, C„ still 
deviates strongly from the expected cubic dependence, 
exhibiting a hump in Cy/T^ which is directly related to 
the so-called boson peak observed by neutron or Raman 
vibrational spectroscopies^. To explain these, it was con- 
sidered that atoms, or groups of atoms, are tunnelling 
between two equilibrium positions, the two minima of a 
double well potential. The model is known as the two 
level system (TLS)*'^. In the standard TLS model, these 
tunnelling excitations are considered as independent, and 
some specific assumptions are made regarding the param- 
eters that characterize them (for a review see for example 
Ref.E). 

In recent years an intriguing magnetic-field dependence 
of the dielectric and coherent properties of some insulat- 
ing glasses was reported. In 1998 Strehlow et al. observed 
a sharp kink at T = 5.84 mK in the dielectric constant 
of the multi-component glass BaO — AI2O3 — Si02^ mea- 
sured in very weak magnetic field, of the order of 10 /iT. 
The effect was several orders of magnitudes larger than 
what is expected, considering the absence of magnetic 
impurities in this insulator. Later studies carried on in 
magnetic field ranging up to 25 T and temperatures well 



below 100 mK revealed for several materials a complex 
and strong dependence of the dielectric response on the 
external magnetic field, on the applied voltage and on 
temperature^i^i^i^. Even more surprising phenomena 
were observed in the spontaneous polarization echoes ex- 
periments on BaO — AI2O3 — Si02^^. The increase of the 
echo amplitude by a factor of 4 was reported when vary- 
ing the magnetic field in the range of to 200 rriT. Sim- 
ilar results were later reported in the case of amorphous 
mixed crystal KBri-xCN^^^. Until now a unified the- 
ory does not exist while some contradictions are present 
in the recent works. We have clear evidence that such in- 
triguing magnetic properties are in connection with the 
low-energy tunnelling excitations present in almost all 
amorphous solids. These tunnelling states are known to 
be proved very well by the echo experiments^. On the 
other hand these excitations have been reported as show- 
ing up at very low temperatures (usually T < 100 mk), 
where the TLS are responsible for the thermal and dy- 
namical properties of glasse s ^^1^^ . 

Several generalizations of the standard TLS model 
have been reported after the anomalous behavior of 
glasses in a magnetic field. The main question for such 
models is how a TLS should interact with the magnetic 
field, being not clear how the tunnelling entity would 
acquire a finite magnetic moment. According to the pro- 
posed solutions, the models can be divided into "orbital" 
and "spin" models. 

In the orbital models, the tunnelling entities are not 
simple two-state systems, but performing some kind of 
circular motions. Due to the presence of the magnetic 
field, a charged particle moves on a loop encloses a 
magnetic flux and thus can acquire an Aharonov-Bohm 
phase. In order to obtain such a closed trajectory, a 
"Mexican- hat potential" was proposed by Kettemann, in- 
stead of the usual double- well potential^^. The resulting 
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flux (j) = nr'^B, with r being the hat radius, proves to 
be by several orders of magnitudes smaller than the flux 
(j)Q = h/e. Even though the existence of cluster configu- 
rations of atoms or molecules containing up to = 200 
units which contribute in a collective tunnelling were re- 
ported, (see for example^^ and references therein) it is 
improbable that such an effect can be extended to an 
amount of the order of TV ~ 10^. A detailed discussion of 
this aspect has been presented in Reflrzl. The combina- 
tion of the Aharonov-Bohm effect with a flip-ffop config- 
uration of the two interacting TLS has been proposed to 
explain the dependence of the echo amplitudes on the ap- 
plied magnetic fieldi^. A different modality to consider 
the occurrence of an apparent flux phase having the cor- 
rect order of magnitude was proposed by Wiirger^^. The 
mechanism proposed consist of pair of coupled TLS and a 
non-linear coupling to the external voltage. Still, the re- 
quired degeneracy of the nearby TLS is not in accordance 
with the known distribution for the tunnel splitting. An- 
other possibility for the formation of closed loops was 
proposed by Le Cocheo^, namely a jagged and uneven 
potential landscape between the two wells. 

To conclude, the "orbital" models can provide an ex- 
planation for some of the magnetic field effects by con- 
sidering the fiux depcndancc of the tunnelling splitting. 
Unfortunately, some assumptions have been made which 
cannot be reconciled with the standard features of the 
tunnelling model. 

The spin modela^S provide an alternative mechanism 
for the observed magnetic properties, especially for the 
polarization echoes^!. A direct coupling between the nu- 
clear spin of the tunnelling entities and the applied mag- 
netic field is considered. We can easily notice that the 
multi-component glasses used in the echo experiments 
contain one or several kind of atoms that carry a nuclear 
quadrupole. For example in the case of BaO — AI2O3 — 
Si02 we find the abundant isotopes Al{I — 5/2) and 
less frequent ^^'^Ba{I = 3/2); for the boro-silicate we 
find the abundant isotope ^^B{I = 3/2). However, the 
echo experiments convincingly put into evidence the role 
of the nuclear quadrupole moments in glasses. Accord- 
ing to this model the magnetic properties should not 
be measured in materials whose nuclei contains no fi- 
nite quadrupole moment. Until now no counter-example 
has been reported. Moreover, recently the role of the 
quadrupole moments have been systematically studied 
and confirmed by the echo experiments in ordinary glyc- 
erol C^HsOz and deuterated one CsDgO^. 

The question, wether or not the nuclear quadrupole 
should be responsible for the magnetic field dependence 
of the dielectric constant seems to be natural. Relevant 
experiments involve temperature where the thermal en- 
ergy is significantly larger than the quadrupole splitting. 
For thermal tunnelling systems, with TLS splitting (E) 
much larger than the quadrupole splitting, it was proven 
already that the resonant, or van Vleck susceptibility is 
rather insensitive to the nuclear quadrupole motion^^. 
For such systems the relaxation or Debye part of the 



susceptibility shows a large magnetic field dependence's. 
By considering TLS with small splitting comparable to 
the nuclear quadrupole energies and using a numerical 
estimation, a pronounced dependence of the electric per- 
mittivity of the applied magnetic field was obtainec^. 

In this article we have studied the magnetic field de- 
pendence of the dielectric susceptibility and the heat ca- 
pacity in cold glasses taking into account the quadrupole 
effects. The general expressions in the mentioned cases 
are complex, however, we have obtained these thermody- 
namic behaviour in the small and large magnetic fields 
limit. In Section II we introduce the nuclear spins in the 
frame of the TLS. The resulting magnetic field dependent 
part of the susceptibility is given in the Section III, us- 
ing the perturbation expansion for two different regimes, 
namely small and large magnetic fields. We end this sec- 
tion with a discussion of our results and a comparison 
with available experimental data and previous calcula- 
tions. In section IV, the specific heat has been calcu- 
lated in three separate terms, namely TLS part; angular 
momentum and nuclear spin parts. Finally the conclu- 
sion is presented with a discussion of our results and a 
comparison with available experimental data3^. 



II. TLS WITH A NUCLEAR SPIN 

The standard TLS can be described as a particle or a 
small group of particles moving in an effective double- well 
potential (DWP). At very low temperatures, only the 
ground states corresponding to the two wells are relevant. 
Using a pseudo-spin representation the Hamiltonian of 
such a TLS is written 

Htls = iAocr^ + ^ACT^, (1) 

where az is the reduced two-state coordinate, 

a, =\L){L\-\R){R\. 

The eigenvalues of cr^ are ±1, labelling the localized 
states in the two wells (left well L, right well R), while 
the tunnelling matrix is taken into account by 

a, =\L){R\ + \R){L\. 

We denoted by A the energy off-set at the bottom of the 
wells, and by Aq the tunnel matrix element. According 
to the randomness of the glassy structure, the energy 
difference between the two wells 

E = ^A2 + A2, (2) 

have a broad distribution. The energy off-set and the 
tunnelling matrix element obey a distribution law 

nA,Ao) = ^, (3) 
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where Pq is a constant. It is useful to define new spin 
operators using the relations 

(Jx = UGE - WUA, 

Gz = WOE + UCTA, (4) 

such that the TLS Hamiltonian becomes diagonal 



D 

TLS 



E 

2 a^. 



We used the notations 



- ^ 
E 

= 1. 



w 



A 

'e 



(5) 



(6) 



which satisfy v? H 

For the moment there is no rigorous theory for tun- 
nelling in glasses. It is assumed that atoms or groups of 
atoms participate in one TLS. As we mentioned before, 
in the case of the multi-component glasses, one or several 
of the tunnelling atoms carry a nuclear magnetic dipole 
and an electric quadrupole. When the system moves from 
one well to another, the atoms change their positions by 
a fraction of an Angstrom. 

We can describe the internal motion of the nuclei by a 
nuclear spin I of absolute value = h^I{I + 1). This is 
related to a magnetic dipole moment gfij^I/h, where g is 
the Lande factor and iin ~ 5 x 

10^27 j/y jg ^j^g nuclear 
magneton. The magnetic dipole couples to an external 
magnetic field B — Bcb and give rise to a Zeeman term 



H 



(1) 



(7) 



where the frequency ez = gfJ-NB/h is directly propor- 
tional to B. 

For a nucleus with spin quantum number / > 1 the 
charge distribution p(r) is not isotropic. Besides the 
charge monopole, an electric quadrupole moment can be 
defined with respect to an axis e 



[3(r-e)2 



p{r). 



(8) 



This can couple to an electric field gradient (EFG) at the 
nuclear position, expressed by the curvature of the crystal 
field potential. The potential describing this coupling 
reads as^l 



Vq = 



-eQ 



1(21-1) 



(9) 



The bases used here (61,62,63) are the principal axes 
of the tensor Vij which describes the electric field gra- 
dient, and 6 is the electron charge. According to the 
Laplace equation the potential obey V11 + V22 + V33 = 0. 
If we define the asymmetry parameter 77 = ^22-^33 ^ ^j^g 
quadrupole potential can be expressed as: 



Vii 



Vq = e,[3lf + r,iL 



I'] 



(10) 



where we denote by Cq 
pling constant. 



ii(2i^-i) ^^"^ quadrupole cou- 



III. STATIC DIELECTRIC SUSCEPTIBILITY 

Our purpose is to determine the magnetic field depen- 
dent part of the static dielectric susceptibility for a TLS 
coupled with a nuclear quadrupole. Therefore we have to 
consider the interaction of the dipole operator {l/2)Paz 
with an external electric field F. The dipole moment 
arises from the relative motion of partial charges related 
to the atoms forming the tunnel system. This interaction 
will be described by 



If we denote 



Vf = {P- F)az 



Z = Tr{6-'^^}, 



(11) 



(12) 



the partition function of the system, where H is the total 
Hamiltonian; the trace, denoted later by < ... > involves 
two-state and spin variables. Here f3 = -jt^, Kb is the 



Boltzmann constant and T is temperature. We can ex- 
press the susceptibility as 



X 



dF^ /' 



(13) 



where / = —^InZ represents the free energy. The sta- 
tistical average implies also an integration over the two 
parameters of the TLS, the energy off-set and the tun- 
nelling matrix element, according to the distribution law 
given by Eq. ([3]), as well as an integration over the nu- 
clear quadrupole parameters; this step, denoted by (...), 
will be discussed at the end of the calculation. We can 
write the partition function in terms of the energy levels 
Ei^m (eigenvalues of the total Hamiltonian) 



(14) 



i=L/B.m=-I 

Therefore, the susceptibility contains three terms: 



X« = 



z ^ 



Qp2 



Z2 



n 2 



dF 



E,. 



E' 



dF 



:E,,_ 



(15) 



The eigen energies {Ei^m) of the full Hamiltonian can not 
be obtained analytically in the general case because the 
Zeeman and nuclear quadrupole terms do not commute. 
In this respect we will consider two different regimes 
where we can obtain the analytic results in a pertur- 
bation frame work. In order to develop a perturbation 
expansion we will consider two separate small 
external magnetic field and a large one, respectively. 
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A. Small magnetic field 

We begin with the case in which the Zeeman term is 
smaller than the nuclear quadrupole potential, ^ Eg, 
thus the Zeeman term is treated as a weak perturbation. 
If we consider the simplest symmetric case, 77 = 0, then 
the nuclear quadrupole potential in the Ii basis has the 
following representation 



\I,m) (/, m|. 



(16) 



where we have defined V/ , 



3m 



/(/ + !), and Eq^^ 



as the quadrupole coupling constant in the left and the 
right well. Therefore, in the (Te basis we can write 

Hq = [e,V/,„i + e'gVi^,n{waE + u(ta)] \I,m){I,m\. (17) 

Here we have defined Sq ~ ^{e^ + e^) and = 
^{Sq — £^)- We can split the total Hamiltonian in two 
parts, 



H = + 5H, 
where the unperturbed part is 



-CJe + EqVl^m - mSz COS f 



(18) 



|/,m)(/,m|,(19) 



and the perturbation term is 

5H = [{e'^Vi,rr+ PF)[waE +u<jA)\\I,m){I,m\ 

- a^' e^\I,m){I,m + l\ 

- af £,|/,m)(/,m-l|. (20) 
We have defined 

^Jl{I + 1) - m{m ± 1) sin O'e'^^^'^''^ 



(21) 



where 9' and (j)' are the direction of magnetic field in the 
basis (ei, 62, 63); here we have assumed that ei = e^. 

We apply the perturbation theory to get the second 
order correction for the energy levels of the system, Ei^m- 
Therefore, the first order correction is 



se: 



(1) 



w{P-F + e'Vi,„^){i\aE\i). 



(22) 



For the sake of simplicity we use {i\(7E\i) instead of 
{i,m\aE\i,m), which is independent of the magnitude of 
m. The second order correction of the energy reads 



6E^ 



(2) 



\{i,m\6H\j,n)\'' 



E. 



(0) 



E 



(0) 



a 



'm\2 



Jm\2 



{2m — \)eq — EzCOsB' {2m + \)eq — SzCosl 



+ -{p.F + e'^Vi,^Y{i\aE\i), (23) 



where Ef'^ is the eigenvalue of the unperturbed part of 
the Hamiltonian 7?^°^. Taking into account the second 
order correction of energy, the susceptibility is calculated 
by Eq. (fT5|) . By using the following relations 



d_ 

OF 



Ei,m\F=o = wP{i\aE\i), 



and 



9F2 



E, 



i,m\F=0 



2PhP 
E 



'i\<TE\i), 



(24) 



we can immediately express the second term of the sus- 
ceptibility 



X2 = PP^W^ 



The third term is 



X3 



(25) 



(26) 



where we defined t = tanh(/3i?/2). The first term of the 
susceptibility contain the second derivative of the energy 
levels with respect to the electric field. After some simple 
algebra xi, can be written as 



2P^u^ 



Xi 



E 



ZE 

m 

where Z = J2,n [e'^'^-'-"^ + e"^^" 
A,„ = £g V/,„i - £zm cos 9' + 



e-^^--] , (27) 
Defining A™ as 



'm|2 



(2m — l)£q — Ez cos9' (2m + l)^^ — Ez cos 9' 



one can write that 



{i\(TE\i) + A„ 



(28) 



The partition function can be written as Z 



E„.^Osh[!^ + ^{E'^Vj,,n) ^ 

term of the susceptibility becomes 



Then the first 



2p2y2 



(29) 

If at this point we assume e'^ = and B = we will 
simply recover the result of standard TLS model, xi = 



1. Thermal expansion 

The experiments we are addressing are performed at 
temperatures between few tens and few hundreds of mil- 
liKelvin, where the thermal energy ksT ~ 10~24j jg 
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few orders of magnitude larger than the Zeeman term 
(at B = IT for example iijsfB/kB ~ OAmK) or the 
quadrupole coupling (echo experiments suggested that 



E 



(l-tW.?x 



,,'m 12 



(2m — l)eq — Ez cos 9' {2m + l)eq — Ez cos 6' 



(30) 



In order to complete the calculation of the magnetic 
field dependent part of the static susceptibility we need 
to perform the integral over the TLS parameters, con- 
sidering the distribution function defined in Eq. Q, for 
this regime (e^ <g; e^), which can be rewritten as 



where k — 



Axi{B) - 

21+1 (Tlir 



2P^K^ 



{2I+ir " e. 



•3 /2 ' 



(31) 



V? 

I"- 



■], and ^ is a constant 



CUm = J2m{{2m-1) ~ (2m+l) 

(nonmagnetic-dependent) which comes from averaging 
over the TLS parameters which depends on the mag- 
nitude of the lower and upper bound of the two-level 
splitting, Aomin and Emax- It can be easily shown that 
in the case of same EFG in the two wells (e^ = e^), the 
dependence of the static dielectric susceptibility on the 
magnetic field will vanish. 



B. Large magnetic field 

Let us denote by (61,62,63) the basis in the left well 
and by (6^, 62, 63) the corresponding one in the right well. 
If we suppose that 63 || 63 || z and e'^, 61 make an angle 
with each other—, we can write the total Hamiltonian in 
the basis of the operator — Iz in the following form 



i?^ + (f + PF) ■ 1 



H^-i^ + PF)-! 



(32) 



77^ are defined for particles We denote £™(f = {m\5V^'^]n) and 



where 1 is the unit matrix of rank {21 -f 1). — Vq 
ff^')^ and ff« = 

in the left and right wells, respectively. The state in 
each well (L or R) can be characterized by | V'/l^'*) =1 
L{R)) I /, to) and consequently the matrix element of 
the Hamiltonian is defined by the following equation. 



After some simple algebra (see for example Ref. [281 ) we 
can easily show that for i?^ „ we obtain 



tz at 200toT, so smaller than one niK). We can H 
conclude that the parameters j3tz, P^q are small, which 
justify an expansion of xi and the exponential factors like 
g±/3ez,g^ Carrying out such an expansion in the thermal 
energy, we can easily show that 



^ = -62 cos 6* m5,n,n+ f-^5m,vym,r. 



tz sm 



^ ^m — 2,n J- m,m — 26 -r ^ ^in+2,n m,m+26 



(34) 

where d and define the direction of magnetic field in 
the basis (61,62,63); a™ and a™ are denoting 



\//(/+l)-?7i(m±l). 



(35) 



and 
T 

T 

^ m,m+2 



1/(/+1)-^to2 + |[/(/ + 1)-3to^] 



1 r 

4 



(3 - ^)\I{I + 1) - (m - l)(m - 2)]5 



X [/(/ + 1) -TO(m- 1)]5 

^ ^(3 - ^)\I{1 + 1) - (m + l)(m + 2)]3 



X [/(/ + 1) - to(to + 1)]2 



(36) 



In a similar manner we can obtain the matrix element 
of the Hamiltonian for the particles in the left well iJ^ „ 
by taking = and changing 1 — > in Eq. (l34|) . 
We will use the perturbation approach assuming that 
2> £q (here for satisfying the perturbation condition 
we consider the special magnetic direction in a way that 
Q ^ f )■ It proves to be useful to change the repre- 
sentation of the Hamiltonian into the basis of the spin 
operators oe and ga- In this basis we split the total 
Hamiltonian of system as B^ -I- SV , where the un- 
perturbed part, will now contain the diagonal matrix ele- 
ments of the total Hamiltonian except those coming from 
Vp; therefore the perturbation, SV ^ will contain the non- 
diagonal terms of the Hamiltonian together with Vp- We 
find that 



m— — I 



E 

<7e — {mtz cos9")\ 



/, to) (/, to| 



and 



{m\5V\n) = Eynn+VmnWaE +VninU(JA- (37) 



£mri — „ (^mn "I" ^rnri) J 



-'mn 



-mn) ■ (38) 



B 



L{R) _ 



^^fm I I <lf ^ (33) 



The diagonal part of the perturbation term is 

Vf = {was +uaA)P ■ F. (39) 
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The first order correction of the energy level will be ob- 
tained as 

^-S'l^m = wP-F(j|cr£;|i)+emm + ?7mmw(j|cr£;|i), (40) 
and the second order one will be 

"^i,m ^ " (0) _ p(0) 

j^i i.ra j-.m 



{JWA\^)\' 



p(0) _ p(0) 
i.m j.m 



c.c. 



E\Vmm\ U -^-y 
p(0) _ p(0) + 2^ 



\{'i-WE\i)VmnW + 

p(0) _ En(O) 
i.m i.n 



(41) 



Since we have assumed that 5V do not contain any diag- 
onal term, emm = ?/mm = 0, the corrections to the energy 
levels simplify to 



^41 = wP-F{i\aE\i) 



(42) 



and 



5E) 



(2) 



E 



{i\(TE\i) 



E 



\{i\(7E\i)'nr,inW + £mnl| 
i.rn i;n 



(43) 



For the large magnetic field regime we calculate the 
first and the second derivatives of with respect to 
the electric field 

d 



and 



^2 



,E, 



i,m\F=0 



i\(JE\i)- 



We can easily express the susceptibility 



ZE 



(44) 



(45) 



Here Z = Xlmle^^^" '" + e~^^"-"'] is the partition func- 
tion, and Ed,m and Eu^m are defined by 

E n 

+ E ;m — ~ ■ (^^^ 



n^ra 



p^"' „ py 



out a thermal expansion in terms of (3eq, (3s z and find 



[1 -I- mcos e"(3ez -pY. 



\{iWE\i)-nninW -f '^mn ^ I 

— COS 6" ez{m — n) 



+0{P^e^), (47) 
After some algebra, we can easily write 



X 



(1) 



-t + 



-(1-Ox 



E E{2I + r 
^ cos9"ez(m — n) 



+ 0(/3V).(48) 



Similar to the small magnetic field regime, from 
Eq. (fT5|) and Eq. (|44| we can straightforwardly show that 



(49) 



The magnetic field dependent correction of the suscepti- 
bility will be 



S^fB) = ^^^^ (l-t^)wu^ ^ 

' E{2i + iy ' 



^ COS 0"ez(m — n) 
(50) 



If we set Eg = 0, we can assume that / = 0, we recover 
once more the result of standard TLS model, 

2P^u^ 



Xtls 



E 



-t + P^w^(3{l-t^). (51) 



By averaging over the TLS parameters according to 
the distribution function ([3]), we obtain: 



5x{B) 



Re[r4„£n 



(2/-fl) ^ cos6"ez(m - n) 



(52) 



where ^' is the numeric constant (nonmagnetic- 
dependent) which depends on the magnitude of E^ax 

and Aomm- We can easily see that 6x{B) ~ where 
£| (X [(4)2 - (e«)2]. 

As we expected from the beginning if we neglect the 
phase difference between the nuclear moment in the wells 
and if we assume that the EFG in the both wells are 
the same (e^ = e^), then we do not find any magnetic 
dependence in the second order perturbation (Eq.[5T]and 
Eq. 



1. Thermal expansion 

Let us remember here that the thermal energy exceeds 
few orders of magnitude the Zeeman energy and the nu- 
clear quadrupole potential. For this reason we can carry 



C. Discussion on the field dependence of dielectric 
susceptibility 

In this section we have addressed the question whether 
or not the coupling of the two-state coordinate and the 
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nuclear spin variables through the nuclear quadrupole po- 
tential in an inhomogeneous crystal field can been taken 
into account as the source of the magnetic field depen- 
dence of the dielectric susceptibility. 

Analyzing the existent data regarding the dielectric 
susceptibility as a function of the magnetic field of differ- 
ent oxide glasses and mixed crystals^i^iSiiiS, we can no- 
tice along with a pronounced bump around 200mT, some 
irregular oscillations. The curvature of the susceptibility 
changes its sign few times up to 5T of the magnetic field. 
The amplitude of the real part of the susceptibility seems 
to be about 10 percent of the TLS susceptibility, corre- 
sponding to approximately 10~^ of the dielectric con- 
stant. Relevant experiments involve temperatures from 
tens of milliKelvin to a few hundreds. We can observe 
that in this range, the dielectric constant vary with the 
inverse temperature, 1/T. 

The present work reconsider the static dielectric sus- 
ceptibility of a glassy systems in an external magnetic 
field implementing a perturbation approach for the en- 
ergy levels . Starting from the coupling of the nuclear 
quadrupole with the tunnelling system, calculations have 
been performed for two different regimes, i.e. small and 
large external magnetic fields with respect to the nuclear 
quadrupole potential. We have found that the magnetic 
field dependent part of the susceptibility in both regimes 
is as the following 

As it is obvious from the above results the dielectric sus- 
ceptibility depends on the magnetic field in the second or- 
der correction of the perturbation scheme via the Zeeman 
and the nuclear quadrupole terms. This correction is di- 
rectly the result of different electric field gradients (EFG) 
in the two wells. In the small magnetic field regime the 
correction increases with while it decreases inversely 
with field (1/-B) in the large magnetic field. The depen- 
dence on the temperature was put into evidence in agree- 
ment with the existent experimental data. The magnetic 
field dependence will disappear in the second order cor- 
rections if we take the same EFG in both wells. We 
would like to remind that the previous existing pertur- 
bation calculations^^ which has been considered the same 
EFG provides a magnetic field dependence of the dielec- 
tric susceptibility only on the fourth order of expansion. 

The theoretical approach we have used do not allow 
us to calculate the static dielectric susceptibility for the 
case where the applied magnetic field has such values 
so that the Zeeman energy and the nuclear quadrupole 
have similar magnitudes, ez ~ e<j. But for this regime 
we know thatSS the magnetic field dependent part of the 
static susceptibility x is proportional to the small ratio 
[{£z£q)/T^]'^ , this ratio being of the order of 10~^ (one 
million part) for T IK. 



We have to mention here that our calculation provides 
only the static limit (w = 0) of this contribution. 

The dynamical susceptibility of a simple TLS was ex- 
pressed by Jackle^^ 

X¥U^) ^ (P^w^ 2!^ , , y (53) 

^^^^ ^ \ cosh^{f3E/2)Lj + ij/ ' ^ ' 

where 7 is the relaxation rate. Now, by considering the 
coupling of the nuclear quadrupole to the tunnelling mo- 
tion, the dynamic susceptibility proves to be more com- 
plicated, but it can be generally expressed asSi 

x'-^'ico) = (p'w' / y-A,„^^\(54) 

' \ cosh2(/3i?/2) "a. + *r™/^ ' 

where represents the eigenvalues of the relaxation 
rate matrix, and Am the corresponding amplitudes (m — 
0, . . . , 2/ -(- 1). In the limit of w — > 0, an approximation 
of the susceptibility is provided by Eqs. ||5T|) and ([5^ for 
the two considered regimes. Reminding that the relevant 
experiments are performed on frequencies in the range 
of kilohertz, and that the relaxation rates should be ex- 
tremely small quantities (details will be give elsewhereSi) . 
We can conclude that such a limiting case cannot provide 
a proper explanation for the existent experimental data. 

The method we developed cannot provide a plausible 
amount for the relative magnitude of the magnetic field 
dependence part with respect to the simple TLS suscep- 
tibility. We do not find oscillations of the dielectric sus- 
ceptibility at magnetic fields of the order of IT. 

IV. THE HEAT CAPACITY 

Our purpose is to find the magnetic field dependence of 
the heat capacity in multicomponent glasses. To the best 
of our knowledge there are not many experimental data 
for heat capacity of (nonmagnetic) glasses in the presence 
of magnetic field. Therefore, we assume that the tunnel- 
ing particles have an angular momentum, J, (intrinsic an- 
gular momentum such as electron spin). Then, we must 
add an additional term to the Hamiltonian which in the 
presence of magnetic field leads to 

H^z'^ = -e'^Jz. (55) 

The energy scale is e'^ = gsfJ'BB/h, fiB is the Bohr mag- 
neton and gE is the electronic Lande factor. Here we 
have neglected the spin-spin interaction in the Hamilto- 
nian. We have also assumed that the quantization axes 
of angular momentum, J^, is in the direction of es (for 
more detail, see the Appendix). The total Hamiltonian 
of the system is then 

H = Htls + Hz + Hq. 

Here Hz = H^z^ + i/^^^ (defined in Eqs. and 

HQ = {^-^)Vii + C-^)Vi. (56) 
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where Vq'"^^ is defined in Eq. (fTO|) for the particles in 
right (left) welP-. The heat capacity is expressed by the 
following relation 



Co, 



1 



KbT^ \ dp 



92 InZ 



(57) 



where Z is the partition function which is a sum over the 
TLS, spin and angular momentum degrees of freedom. 
If we neglect the phase difference between the nuclear 
moments in the two wells (taking — Q) and assume 
that the EFG in both wells are the same (e^ = = e^), 
then the partition function can be written in the following 
form 



Z ^ Z' 



TLS 



SS' 



(58) 



where 



Ztls 
S 

S' 



Trje"''^^^-''} 2cosh(/3S/2), 



Trje-'^^i^'}. 



(59) 



Therefore, according to Eq. ([57|l and Eq. ([58| . we find 
that 



1 



KbT^ \ a/32 



a/32 



a/32 



(60) 



By some calculations and after averaging over TLS pa- 
rameters, we can easily show that [see for example 
Ref. dil)] 



G 



TLS 



6 



(61) 



which is independent of the magnetic field. The Zeeman 
contribution of the tunneling particle's angular momen- 
tum (J) in the partition function is simplified to the fol- 
lowing form, 



S' 



E 



sinh[/3e',(J-t-i)] 



m'=-,/ sinh[^] 



(62) 



Thus the Zeeman contribution of the angular momentum 
in the specific heat is given by 



csch 



P' /d'^lnS' 



KbT^ 

2 I P^'z 



a/32 



2_/2 



{2J+l)^csch^ 



X P' 



/34(2J+1) 



(63) 



where P' is the concentration of intrinsic angular momen- 
tum times a constant which comes from the averaging 
over the TLS parameters. 



Finally, the contribution of the remaining term of the 
Heat Capacity to the specific heat, C^-^\ is given by 



a2 In 5' 

a/32 

T' - H-(l) 



(64) 



where S = Tr{e-P"'} and H' = H^' + Hq. 

As we have mentioned in previous section the ther- 
mal energy exceeds few orders of magnitude the Zeeman 
energy and the nuclear quadrupole potential. There- 
fore, /3ez and [3eq are very small quantities, consequently 
j3H' -C 1. Therefore, Eq. ([M]) is approximated by 



^B/3^^1n 



Tr(l - PH' 



2 



+ ...) 



21+1 



TriH'^) — -) 



(65) 



and finally we arrive at the following expression 

C(^) - P"KBP^7iel+l2e^g]+0{(3^e^). (66) 

where P" is the concentration of the nuclear spin times 
a constant which comes from the averaging over the TLS 
parameters. Here 71 and 72 are numerical constant which 
have been defined by (see appendix A), 



71 



72 



21 



21 



^ [m2 cos2 0' + sin2 0' [/(/ + 1) 

m 

— VfT2 



y2 

m,m — 2 



2/-f 1' 



(67) 



Eq. ([66]) shows a quadratic dependence on the mag- 
netic field. In the large field regime C*^^^^ contributes 
as a constant value to the whole specific heat. Thus, the 
field dependence of the specific heat in the high magnetic 
fields comes from C^-^^ which shows a quadratic behaviour 
versus the magnetic field. However, for small field regime 
we have Ez ^ which makes C^^ to have the dominant 
effect {C^''> < C^^^)). It can be shown that when the 
magnetic field is small the leading term of C^' is given 
by the following expression 



,,J{J+l) 



which shows behaviour for low magnetic fields. 



(68) 



A. Discussion on the field dependence of the Heat 
Capacity 

In this section, we have studied the magnetic field de- 
pendence of the Heat capacity of glasses at low tempera- 
ture. We have assumed that the tunneling particles carry 
both an angular momentum and a nuclear spin. So, the 
magnetic field couple to both terms. 



9 



0.8 




5 10 15 

Pe,/2 

FIG. 1: Magnetic field (thermal) variation of the C' 



We have examined our results numerically for arbi- 
trary angle of EFG and different phases between the 
quadrupole moments on each well. The final results are 
not influenced if we only consider the simple case of = 
and considering that — — eq. Thus, we have re- 
ported our results in the mentioned special case where 
we can find an analytical expression for the specific heat 
behaviour. In this respect, the total partition function 
of system can be splitted into the multiplication of three 
separate terms. Therefore, the Heat capacity is com- 
posed of three terms, the TLS contribution (Ctls), the 
nuclear (C^-^^) and the angular momentum (C^^^^) parts. 
The magnetic field dependence will arise from C^^^ and 

We have introduced the generalized Hamiltonian which 
takes into account the mentioned degrees of freedom. 
The direction of external magnetic (6*', (j>') field has been 
considered arbitrary with respect to the EFG direction 
of each TLS. However, the outcome will be averaged over 
this spherical angle which can be contracted as a constant 
to the final results. We have calculated the contribution 
of Ctls and C'^^^ exactly while the effect of nuclear spin 
(C^^^ ) has been treated in a thermal expansion up to the 
second order corrections. 

The different scale of Zeeman energy for the nuclear 
moment and the angular momentum of the tunneling en- 
tity is roughly Cz/e^ ~ 10^'^. Thus, for low magnetic 
fields /Se'j, < 10 the dominant effect comes from the angu- 
lar momentum part {C'^^^^). Our calculation shows that 
at low magnetic field this contribution is proportional 
to the square of magnetic field, C^'^'^ ~ B'^ fEq. ([S7)) . 
However, C*-^^-* reaches a maximum at fie'^ ~ 3 and then 
decreases to zero for /3e^ > 20 (see Fig.dl])). 

For low magnetic field {B < lOT) our result predict a 
Schottky like peak in the specific heat of glasses which 



have nonzero angular momentum for the tunneling entity. 
Or at least the sample has some impurity with nonzero 
J (angular momentum). This result is in agreement with 
the recent experimental observations for Duran, AlBaSi^ 
and SuprasiP^. The data show explicitly an upward in- 
crease of the specific heat versus the magnetic filed for 
B < 0.5T and monotonic decrease for i? > IT at T < IK. 

We have also predicted that for high magnetic fields 
(/Se'j, > 20) the dominant contribution to the specific heat 
comes from the nuclear moments which shows a quadratic 
dependence on the field (Eg. 

In this paper we have studied the special case that 
the impurities have only the intrinsic angular momentum 
part. For the general case of impurities (for example 
electrons) with angular momentum, J > ^ (spin + orbit), 
we should take into account the effect of spin - orbit 
interaction. Which means that at zero magnetic field the 
levels m' — — J, . . . , J — 1 , J may have different energies 
due to the crystal field splitting. 

However, if we consider J = which means to have 
a pure glass without magnetic impurities like BK7 the 
nuclear quadrupole term is the only dominant term in 
the specific heat. So, the present calculations predict a 
quadratic dependence on the magnetic field, and l/T^ 
dependence on temperature. As far as we know there is 
no experimental result for this case. Therefore it might 
be a good suggestion for future experiments. 



V. SUMMARY 

Based on the polarization echoes experiments— we 
believe that the nuclear quadrupole model should pro- 
vide the rightful explanation for the intriguing behav- 
ior of the dielectric properties of glasses. A relaxation 
spectrum rather sensitive to the orientation of the nu- 
clear quadrupoles is providing a relaxation susceptibility 
strongly dependent on the applied magnetic fieldSi, driv- 
ing us to the conclusion that the observed magnetic prop- 
erties of glasses might have a relaxation origin. Other 
effects like non-linearities with the voltage^, or cooper- 
ative TLS behavior might as well be found contributing. 
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APPENDIX A: GENERALIZED TLS 
HAMILTONIAN WITH BOTH NUCLEAR SPIN 
AND ANGULAR MOMENTUM 

The aim of this appendix is to find the generahzed 
TLS Hamiltonian with both nuclear spin and angular 
momentum. In the similar way leading to the Eq. 1321 we 
can write the total Hamiltonian as: 



H 



^ I A 1 Aq 1 

. 2 2 

Ao -1 

2 



(Al) 



particles in the left and right wells, respectively. The 
state in each well (L or R) can be characterized by 

I ^^!m],[J,m']) ==1 I J,m') and conse- 

quently the matrix element of the Hamiltonian is defined 
by the following equation, 



IjLiR) _ I iMR) I rrL{R) | MR) . , 

■"mm'.nn' " W[I ,m],[J,m'] I I ^[I ,n],[J,n']> ' ^-^^1 



where 1 is the unit matrix of rank (2/ + 1) x (2J + 1) 

= + Hk and ^ VM 



H§ are defined for 



After some simple algebra we can easily show that for 
H^m'nn' ^6 obtain 



mm' ,nn' 



-26 



^m — 2,n ~r ^ m,m+2^ 



-2ie i 



cos 9 m + e^Tr 



5m+2,n\ 



— e, smf 



j0 msj I -10 msj 



(A3) 



where 6 and (p define the direction of magnetic field 
in the basis (61,62,63); aj? and T^.m' are given by 
Eqs. pSI & I36p . In a similar manner we can obtain the 



matrix element of the Hamiltonian for the particles in the 
left well i?^ „ by taking 6 = and changing 1 — > 
in the above equation. 
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